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ABSTRACT 

Gamma-ray burst X-ray flares are believed to mark the late time activity of the central 
engine. We compute the temporal evolution of the average flare luminosity (L) in the 
common rest frame energy band of 44 GRBs taken from the large Swift 5-years data 
base. Our work highlights the importance of a proper consideration of the threshold of 
detection of flares against the contemporaneous continuous X-ray emission. In the time 
interval 30s < t < 1000s we find (L) oc £- 2 '7±0.i. this implies that the flare isotropic 
energy scaling is Ei so , fl arc oc t~ . The decay of the continuum underlying the flare 
emission closely tracks the average flare luminosity evolution, with a typical flare to 
steep-decay luminosity ratio which is Lfl are /i s teep = 4.7: this suggests that flares and 
continuum emission are deeply related to one another. We infer on the progenitor 
properties considering different models. According to the hyper-accreting black hole 
scenario, the average flare luminosity scaling can be obtained in the case of rapid 
accretion (i acc <C t) or when the last ~ 0.5Mq of the original 14Mq progenitor star 
are accreted. Alternatively, the steep oc t~ 2 - 7 behaviour could be triggered by a rapid 
outward expansion of an accretion shock in the material feeding a convective disk. If 
instead we assume the engine to be a rapidly spinning magnetar, then its rotational 
energy can be extracted to power a jet whose luminosity is likely to be between the 
monopole (L oc e~ 2t ) and dipole (L oc i~ 2 ) cases. In both scenarios we suggest the 
variability, which is the main signature of the flaring activity, to be established as a 
consequence of different kinds of instabilities. 
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1 INTRODUCTION 

Observations of long gamma-ray bursts (GRBs) by a number 
of spacecrafts during the last decades revealed that their 7- 
ray prompt emissio n is over after ~ 10— 100 s. However, after 
the launch of Swift l|Gehrels et al .1 12004') evidence is accumu- 
lating that the central engine of the GRB sources could still 
be active hours after the main burst. This is mainly linked 
to the discovery of large amplitude, episodic re-brightenings 
with typical At/t ~ 0.2 in the ~ 33% of X -ray afterglows: 
the X-ray flares (see IChincarmT et al.ll20ld. C10 hereafter, 
for a recent compilation). The temporal properties of X-ray 
flares make it difficult to interpret the observed emission 
in the framework of the external shock scenario (see e.g. 
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lLazzati fc Pernal 2007). M oreover, the strict analogy found 
bv lMargutti et all (|2010bl) (M10 hereafter) between the tem- 
poral and spectral behaviour of X-ray flares and prompt 
emission pulses strongly suggests a common, internal ori- 
gin. The direct implication is that flares directly trace the 
activity of the central engine. 

To account for the details of the flare emission, GRB 
central engines (CE) are required to have the following basic 
properties (C10, M10 and references therein): 

(i) Energetics. The average flare fluence is ~ 10% the 
prompt 15-150 keV fluence. However, CEs should be able 
to produce X-ray flares whose fluence is able to compete 
in some cases (e.g. GRB 050502B) with the prompt 7-ray 
fluence. The typical early-time flare isotropic energy is ~ 
10 51 erg. The flare peak luminosity decays as L pcak oc t~ 2 ' 7 
for t < 1000 s. 
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(ii) Lifetime. CEs are required to be long-lived and in 
particular to be active up to t ~ 10 4 — 10 5 s after the end 
of the prompt emission. The probability of a CE to re-start 
decays with time. 

(iii) Variability. CEs are required to store and release en- 
ergy in the form of erratic, short-lived episodes of emission 
giving rise to flares whose duration linearly grows with time. 

Ideas on how to revive the CE have been explored 
by a number of authors: flares could be the result of a 
two-stage stellar colla pse leading to cor e fragmentation and 
subsequent accretion l|King et al.l 120051 ). Alternatively, the 
fragmentation of the outer part of an hyper-accreting disk 
into rings of material as a result of gravitational instabil- 
ities could potenti ally lead to large-amplitude va riations 
of th e CE output (IPerna. Armitage fc Zhand 120061 ). How- 
ever, iPiro fc Pfha3 \200H ) noticed that tidal disruption of a 
fragment and accretion by the central black hole might be 
too rapid to ac c ount for the durations of observed flares. 
IProga fc Zhand (|2006l ) proposed that magnetic flux accu- 
mulated around the accretor could stop and then restart the 
CE re lease of energy possibly giving rise to flares. iDai et ail 
(2006) invoked the presence of a differentially rotating mil- 
lisecond pulsar where magnetic reconnection events are re- 
sponsible for the observe d flare emi s sion. M agnetic reconnec- 
tion is also suggested bv lGianniosI (2006) in the framework 
of the internal plus external shock scenario with strongly 
magnetised ejecta: in this case, no late-time CE activity is 
needed and flares come from the revival of MHD instabilities 
triggered by the deceleration of the magnetised ejecta in the 
external medium. 

The entire list of flare models basically refers to two dif- 
ferent C Es: an accretion disk around a newly formed black 
Woos lcv 1999) or a rapidly rotating 



Table 1. Sample of 44 GRBs observed in the period April 2005 
- February 2010 by Sim/i-XRT with X-ray flaring activity and 
measured redshift. 



hole (e.g. iMacFadven 



magnetar (e.g. Usovl 19921: iThompsor] 1 19941 : IWheeler et alj 



2000; iThompson et al.ll2004l ~ In magnetar models the ulti- 



mate source of energy is represented by the magnetar spin 
and flares are likely to be connected to episodes of mag- 
netic energy dissipation. In the case of accretion models, 
the source of energy powering the flare emission is possibly 
of gravitational origin. However, magnetic effects could still 
play the key role, deeply affecting the accretion dynamics 
and release of energy (e.g. iFan et al.ll2005l : IProga fc Zhangl 
l2006l : lLei et al.ll2009l ; lzhang fc Daill2009l ). X-rav flares mark 
the late time activity of the CE after the main power emis- 
sion and provide information about the way the CE power 
is progressively declining. In the following we show that this 
information turns out to be of fundamental importance in 
constraining existing CE models (accretion or magnetar). 

The main goal of this work is to compute the evolu- 
tion of the average luminosity of the flaring component with 
time, taking advantage from the large Swift 5-years data- 
base. A variety of physical mechanisms which might prop- 
erly account for the new episodes of energy release are dis- 
cussed in the context of both accretion and magnetar CE 
models. A pr e vious attempt in this direction was made by 
lLazzati ~et al. (2008) (L08 hereafter): from a sample of 9 long 
GRBs and 1 short GRB (which was the largest sample with 
known redshift available at that time) these authors showed 
that the average flare luminosity follows a power-law decay 
(L)h08 <x t" 1 ' 5 * - 2 and concluded that accretion onto a com- 
pact object could reasonably account for the (L) decay. In 
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the present work we aim at presenting a substantial update 
of the L08 sample, including 44 GRBs with redshift; the 
flare luminosity of each GRB is computed in the rest frame 
energy band which is common to the entire sample: this 
assures (L) to be obtained from strictly comparable proper- 
ties of each GRB. Moreover, no assumption is made on the 
functional form of the flare temporal profile. Finally, since 
flares are detected as sharp features superimposed to the 
X-ray afterglow, it is of fundamental importance to under- 
stand the associated observational biasses: we compute the 
temporal evolution of the flare detection threshold for each 
GRB, and quantitatively discuss the influence of the flare 
detection threshold on the obtained (L) temporal evolution. 

This work is organised as follows: the sample selection 
and data reduction is detailed in Sect. [2] while the compu- 
tation of the flare and underlying continuum components of 
each event is presented in Sect. [3] Section [4] describes the 
method which leads to the definition of a flare detection 
threshold of each GRB. A luminosity-time plane avoidance 
region is presented in Sect. [5] Results are discussed in Sect. 

possible physical interpretations in the context of both 
accretion and magnetar models of the central GRB engine 
are critically analysed in Sect. 16. Ol and lF. 1.21 Constraints to 
physical scenarios able to account for the erratic behaviour 
of the flare emission are derived in Sect. 16.21 Conclusions are 
drawn in Sec. [7] 

The phenomenology of the bursts is presented in the 
rest frame unless otherwise stated. Uncertainties are quoted 
at the 68% confidence level (c.L): a warning is added if it 
is not the case. Standard cosmological quantities have been 
adopted: H = 70Kms~ 1 Mpc- 1 , n A = 0.7, Q M = 0.3. 



2 SAMPLE SELECTION AND DATA 
REDUCTION 

We select the GRBs which visually show flaring activity 
superimposed to the smooth X-ray afterglow, good data 
qualit y and redshift, detect ed by the Swift X-ray telescope 
(XRT, H urrows et al.ll2005h from the beginning of the mis- 
sion up to the end of February, 2010. The sample comprises 
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44 GRBs (TableHJ), with redshift ranging from 0.3 to 6.33. 
XRT data have been processed with the HEASOFT package 
v. 6.6.1 and corresponding calibration files: standard filter- 
ing and screening criteria have been applied. Piled-up Win- 
dow Timing (W T) data have been c orrected following the 
prescriptions bv lRomano et al.l (|2006l ). while piled-up Pho- 
ton Counting (PC) data have been extracted from an annu- 
lar region whose inner radius has been derived comparing 
the observed to the nominal point spre ad function (PSF, 
iMoretti et ai1l2005l ; IVaughan et ail 120061 ). The background 
is estimated from a source-free portion of the sky and then 
subtracted. The 0.3-10 keV background subtracted, PSF and 
vignetting corrected light-curve of each GRB has been re- 
binned so as to assure a minimum signal-to-noise (SN) equal 
to 4. The count-rate light-curves are calibrated into luminos- 
ity light-curves using a ti me dependent count-to-fl ux conver- 
sion factor as described in lMargutti et al.l (|2010al ). This pro- 
cedure produces luminosity curves where the possible spec- 
tral evolution of the source is properly taken into account. 
Each 0.3-10 keV XRT light-curve is calibrated in the com- 
mon rest frame energy band defined by the redshift distribu- 
tion of the sample which turns out to be 2.2-14.4 keV. This 
allows us to make a direct comparison between the X-ray 
afterglows of different bursts while avoiding extrapolation 
of the signal to an unobserved energy banq 2 ]. 
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Figure 1. X-ray afterglow of GRB 060607A in the rest frame en- 
ergy band 2.2-14.4 keV. Red dot-dashed line: best estimate of the 
continuum underlying the flaring emission obtained as detailed 
in Sec. \3\ Blue solid line: 2a threshold of flare detection (i.e. the 
minimum peak luminosity that a flare must have to be detected 
as statistically significant fluctuation superimposed to the X-ray 
continuum), calculated as described in Sec. [J] In this case, for 
t < 1000 s, the threshold is found to be ~ 50% the continuum. 
Inset: residuals with respect to the continuum. 



3 CONTINUUM AND EXCESS ESTIMATION 

We produced a software to automatically identify the 
smooth continuum underlying the X-ray afterglow of GRBs 
with superimposed flaring activity. The procedure is based 
on the x 2 statistics and can be described as a two-step pro- 
cess: 

(i) Identify the continuum afterglow component underly- 
ing the flaring activity; 

(ii) Calculate the flaring contribution as excess with re- 
spect to the continuum. 

A first blind fit of the entire X-ray afterglow (continuum plus 
flares) is done in log-log units using a power-law, smoothly 
joint broken power-law or double broken power-law models. 
If the P-value (probability of obtaining a result at least as 
extreme as the one that is actually observed) associated to 
this fit is lower than 5%, the data point with the largest pos- 
itive residual is removed from the light-curve and a new fit is 
performed. This process is repeated until a P-value > 5% is 
obtained. The F-test is used to choose between the different 
nested models when necessary. The best fitting model satis- 
fying the P-value condition is identified with the underlying 
continuum associated to a particular GRB X-ray afterglow 
(red dot-dashed line in Fig. Q] for the case of GRB 060607A) 



1 GRB 090423 with z=8.3 was not included in spite of the clear 
X-ray flaring activity to have a wider sample common rest frame 
X-ray energy band. 

2 An exception is GRB090417B with z=0.345: this burst was 
added later to the sample and required the extrapolation of the 
observed flux from 2.2-13.5 keV rest frame to 2.2-14.4 keV rest 
frame. 



and is subtracted from the original light-curve. The result- 
ing residuals and respective errors (both the statistical un- 
certainty associated to the original light-curve bins and the 
one coming from the continuum estimation^ are properly 
taken into account and propagated) constitute the candi- 
date X-ray flaring component associated to a GRB: this is 
shown in the inset of Fig. [Qfor GRB 060607A. Note that no 
particular flare functional shape is assumed in this analy- 
sis. This method allows us to account for small variations 
superimposed to the continuum: in C10 flares were instead 
visually identified and then fitted with a specific profile. 

The average GRB 2.2-14.4 keV rest frame flaring com- 
ponent for any time interval ti — tf is computed as: 

(L(t)) 2 ' 2 - 14 - 4koV = 10^ E ^ Lo9(Li) (1) 

where N is the number of GRB displaying a positive flaring 
component at a minimum 2a significance during ti — tf, Li 
is the 2.2-14.4 keV luminosity of the flaring component of 
the i th burst evaluated at at t = (ti +tf)/2. Linear interpo- 
lation is used when necessary. The uncertainty affecting (L) 
is found through standard propagation of the uncertainties 
affecting each flaring component, as determined in the previ- 
ous paragraph. The typical uncertainty on each GRB flaring 
component data point is of the order of 30% the value of the 
GRB afterglow light-curve before the flare subtraction. An 
example is shown in Fig.[T] inset. The use of the logarithmic 
mean in Eq. [T] prevents (L) from being dominated by a sin- 
gle bright event (the use of the median has been proven to 
lead to very similar results); instead, the linear mean of the 
excesses leads to (L) values biassed towards the bright end 
of the luminosity distribution of the excesses at any time t, 
and is therefore discarded. GRB 050904 shows an extended 

3 We take into consideration the complete information coming 
from the covariance matrix of the fit. 
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Figure 2. Panel (a): average flaring component luminosity curve 
of the 44 GRBs computed in the common rest frame 2.2-14.4 keV 
energy band. In the time interval t < 1000 s (t > 1000 s) the best 
fit reads {L/ergs^ 1 ) = IO 54 - 5 * - 1 (t/s)- 2 - 7 * - 1 ((L/ergs- 1 ) = 
10 51,4±01 (t/s)^ 1 - 2 ^ ' 1 ). Panel (b): average re-normalised lumi- 
nosity curve: the luminosity has been re-scaled for the average 
prompt 15-350 keV luminosity observed during the Tqq. The x- 
axis is expressed in T90 units. In the time interval 1 < t* < 100 
the re-normalised flare luminosity follows a power-law decay with 
best fitting model: (-L/^prompt) = IO 1 - 1 * ' 1 (t*)-l-2±0.1 w here 
t*=t/T go . 



flaring activity between 500 s and 3.1 ks: given the high lu- 
minosity of the flaring component and the peculiarity of a 
flaring emission covering the entire plateau phase duration, 
GRB 050904 is not included in the calculation above. This 
guarantees that (L) is not biassed by the contribution of 
a single peculiar burst: however, the inclusion of this event 
would not change the main conclusions of this work. The 
flare luminosity curve of GRB 050904 is investigated and 
discussed in Sec. [6] The result of the application of Eq. [I] is 
shown in Fig. [2] panel (a): in the time interval 15 — 30 s the 
average light-curve of the flaring component is likely to track 
the prompt gamma-ray emission, while for 30 < t < 1000 
s the average light-curve is well represented by a power- 
law with best-fitting index ai = -2.7 ± 0.1. For t > 1000 
s the best-fitting power-law index reads «2 = —1.2 ± 0.1. 
We note that a simple power-law modelling of the entire 
15 s < t < 10 6 s temporal window yields a statistically un- 
acceptable fit (x 2 /dof = 1271.4/50, P-value< 10~ 10 ) that 



systematically overestimates the computed average flare lu- 
minosity for t > 1 ks. The two data points at t ~ 10 3 s 
largely above the detection threshold (green thick area, see 
Sec. [4] for details) are mainly due to the flaring contribu- 
tion of GRB 081028, GRB 090809 and GRB090417B: these 
3 events show a large flare at this rest frame epoch. The 
exclusion of these 3 bursts from Eq. [1] would bring the two 
data points slightly above the threshold. The shown error 
bars account for the la uncertainty affecting the average 
luminosity value: however, the dispersion of the sample lu- 
minosity values at any given time t could be remarkably 
different. For this reason a yellow line- Allied area marking 
the 90% sample dispersion has been added to Fig. [2] both 
panels. 

The GRB prompt 7-ray emission is known to span dif- 
ferent orders of magnitude in luminosity. We correct for the 
different prompt luminosity characterising different events 
re-normalising the flaring component by the average 15- 
350 keV luminosity observed during the BAT T90 (inter- 
val of time of emission of 90% of the BAT fluence). The 
re- normalised flaring luminosity (L*) = (L / L plolrl pt) in the 
time interval t* — t* s is computed as: 



{£•(«*)> 



2.2-14.4 keV 



= 10" 



(2) 



where L* = Li/(Li n ), being (Li, 7 ) the average 15-350 keV 
prompt luminosity during the T90 of the i th burst. For each 
burst, the time t* is expressed in T 90 unitfl Both the BAT 
T90 and the fluence information are retrieved from the burst 
BAT-refined circular^. Figure [5] panel (b ), shows the re- 
sulting re-normalised average flaring luminosity curve: for 
t* < 1 the XRT caught the X-ray counterpart of the prompt 
7-ray emission; in the time interval 1 < t* < 100 the re- 
normalised flaring activity follows a power-law decay with 
best-fitting index a = —1.2 ± 0.1. 



4 FLARE DETECTION THRESHOLD 

Flares are always detected as emission components superim- 
posed on top of the overall afterglow decay: this means that 
whatever the origin of the continuum emission is, the sensi- 
tivity to temporal fluctuations is degraded as the flux level 
-and consequently the statisti cs- of the un derlying contin- 
uum decays with time (see e.g. Morris 2008). This naturally 
requires flares at late times to be longer in duration and/or 
have a higher flux contrast with respect to the continuum in 
order to be detected as statistically significant fluctuations. 
It is therefore of primary importance to understand at which 
level the evolution with time of the average flaring compo- 
nent portrayed in Fig. [2] is a by-product of the decay of the 
underlying X-ray afterglow which acts as a time- variable de- 
tection threshold. 

To this end a set of simulations is run for each GRB of 
our sample: the aim is to determine the flux and luminos- 
ity detection threshold of the flaring component at any time 



4 GRB 051022 was detected by HETE and consequently has no 
measured BAT Tgo: for this reason it is not included in the cal- 
culation. GRB 050904 is also excluded because of its peculiarity 
(see Sect. [OJ - 

5 For GRB 060424 we refer to iRomano et al.l [2006). 



Average GRB X-ray flaring activity 5 




10 8 10 3 10 4 10 5 1 10 100 1000 10000 

Rest Frame Time (s) Time/T 90 



Figure 3. Panel (a): X-ray flaring activity contour plots in the luminosity-time plane. The area is color coded as defined in the middle 
panel. The luminosity is computed in the rest frame 2.2-14.4 keV energy band. The flaring activity of 90% of the GRBs in our sample 
lies below the red solid line. Yellow dashed line: power-law with index -1 plotted for comparison. Panel (b): the luminosity is expressed 
in units of the average 15-350 keV prompt luminosity during the BAT Tgo; times are in Tgo units while the time is expressed in units of 
BATT go . 



t given the observed X-ray afterglow statistics. The con- 
tinuum contribution to each afterglow and associated un- 
certainty have been calculated in Sect. [3] we evaluate the 
continuum at the original light-curve central time bins and 
combine the systematic uncertainty coming from the con- 
tinuum evaluation procedure with a statistical uncertainty 
equal to the relative error on the original luminosity bins of 
that particular burst. A grid of 545 flare peak time values t p k 
equally spaced in logarithmic units between 1 s and 10 7 s is 
generated together with a grid of 171 flare peak luminosity 
values in the interval 10 40 — 10 57 ergs -1 . For each t p k, we 
generate 171 flares with growing peak luminosity values: a 
iNorris et all (|2005l ) profile specified in terms of amplitude A, 
width w and peak time t p k is assumed. The Norris 2005 pro- 
file is defined as / = f(A, Ti,T2,t s ) where t\ and T2 are two 
flare shape parameters related to the rise and decay phases, 
respectively, while t s defines the pulse starting time. We take 
advantage of previous studies on the flare phenomenology 
to reduce the flare parameter space and write n ~ 4uj/3; 
t 2 ~ w/3; t a ~ t pk - 2w/3; w ~ i Q -o.94 t (i-i2±o.i8) ( see C10 
and M10 for details). The signal of the fake flare is then inte- 
grated in the original light-curve bins and a statistical error 
is assigned equal to the relative uncertainty affecting the 
original luminosity curve at the same time. Finally, the flare 
contribution is added to the continuum and the uncertainties 
propagated into a final fake GRB afterglow. To simulate the 
entire procedure, the fake afterglow is then processed by the 
automatic excess estimation routine described in the previ- 
ous section: a new continuum is determined together with 
the associated flaring component, as before. If the newly de- 
termined flaring component contains a positive fluctuation 
with a minimum 2a significance, the peak luminosity of the 
fake flare is recorded as 2a flare detection threshold at t p k 
and the following peak time value is considered. The entire 
procedure is repeated for each t p k , every time starting from 
the minimum peak luminosity value of the grid. The result 
is shown in Fig. [T] for GRB 060607A taken as an example 
(blue solid line): due to the degrading quality of the signal, 



a growing flare-to-continuum flux contrast is required for the 
flare to be detected at later times. 

The 2a detection luminosity curve of the different GRBs 
are then combined to produce the average detection thresh- 
old (Lthit)) in the rest frame band 2.2-14.4 keV. In strict 
analogy with Eq. [T] (L t h{t)} is defined as: 

(L th (t)> 2 ' 2 " 14 - 4kcV = io^E&i^m*) ( 3 ) 

being N the number of GRBs displaying a significant posi- 
tive flaring component in the time period ti — tf, ii,tft is the 
luminosity of the 2a detection threshold of the i th GRB at 
time t — (ti+tf)/2. The re-normalised average 2a threshold 
(L* h (t*)) is computed as: 

(L*(t*)) 2 - 2 - 14 ' 4koV = wfrEtLiLog(Ll th ) ^ 

where L* th = L itt h/{L i>rl ). As before, N collects the GRBs 
with significant excesses in t* —t%. 

The result is portrayed in Fig. [2] in both panels a thick 
green area marks the 2a detection threshold region. From 
panel (a) it is clear that the —1.2 slope is mainly due to the 
selection of bright fluctuations able to overshine the thresh- 
old. The fact that (L(t)) tracks the {L t h(t)) temporal be- 
haviour suggests that the sample of significant excesses de- 
tected in the time period 2-3000 ks are not representative of 
the real population of flares at those times. The bright end 
of the flare luminosity distribution has been sampled and 
a biassed average flare luminosity has been re-constructed. 
Notably, the decay of the detection threshold at t ~ 10 4 s 
allows us to appreciate a steeper decay (a ~ —2.5) of 3a 
significant positive fluctuations extending to t ~ 6 x 10 4 s. 
No more than 2a significant flaring activity can be quoted 
for (1 < t < 6) X 10 5 s due to the flatter decay (and actually 
enhancement) of {Lthit))- The unbiassed average flare lumi- 
nosity curve is therefore likely to be steeper than —1.2 at late 
times. At early times the situation is different: for t < 300 s 
the (L(t)) oc t~ 2 ' 7 decay is not affected by threshold effects; 
the extension of the same temporal law up to t ~ 1000 s 
suggests that (L t h(t)) is likely to play a minor role in the 
time interval (0.3 < t < 1) ks, as well. Finally, Fig. [2] panel 
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(b) shows that the re-normalised average flare luminosity 
tracks the threshold starting from t* ~ 100. (L*(t*)) is con- 
sistent with the —1.2 power-law decay up to t* ~ 800. For 
t* 5s 1000 the re-normalised activity significance is < 2a. 



5 LUMINOSITY-TIME PLANE: THE FLARE 
AVOIDANCE REGION 

The 2a detection threshold luminosity curves calculated in 
Sec. [4] define for each GRB, for each time t, the minimum 
peak luminosity which a flare must have to be detected as 
statistically significant positive fluctuation above the X-ray 
afterglow (see Fig. [TJ blue solid line). We combine this infor- 
mation with the residual luminosity curves derived in Sec. 
[3] to draw an upper limit to the X-ray flaring luminosity as- 
sociated to each afterglow at any time t, as follows: for each 
event we substitute the luminosity value of non significant 
fluctuations with the respective 2a upper limits from the 
corresponding threshold luminosity curve, keeping the sig- 
nificant fluctuations unchanged. In this way we account for 
the presence of undetectable flaring activity. The resulting 
curve divides the luminosity-time plane into two halves: the 
real flaring luminosity curve is expected to lie below this line 
at a minimum 2a c.L: we refer to this region as the X-ray 
flare-area associated to a particular GRB. This procedure 
is repeated for each burst of the sample for both luminosity 
and re-normalised luminosity units. The different flare-areas 
are then superimposed and a luminosity-time plane contour 
plot is created (Fig. [3J: the different colours refer to the 
different number of GRB flare-areas superimposed on a par- 
ticular (ti,Li) pixel. In both panels the flaring luminosity 
of the 90% of the GRBs of our sample lies below the red 
solid lin^E We refer to the area above the red line as the 
luminosity-time plane flare avoidance region. Note that a 
direct implication of Fig. [3] panel (a) is that in the time 
period 60s < t < 400s the flare luminosity function decay 
is steeper than oc a t 90% confidence. 



6 DISCUSSION 

In the previous sections we showed that: 

• For 15 s < t < 30 s the average flare luminosity function 
is likely to track the prompt 7-ray emission; 

• In the time interval 30 s < t < 1000 s we find (L) oc 
t~ 2 No threshold-related argument can be invoked to 
explain this scaling. 

• For t > 1000 s (L) closely follows the temporal be- 
haviour of the threshold of detection: (L) oc i~ 1,2±01 . For 
this reason we believe the un-biassed average flare luminos- 
ity function at this epoch to be steeper. 

Notably, a s i milar L p k — t p k scaling was obtained by 
IChincarini et all (|2010h for t < 1000 s once the intrinsic 
scatter of the relation was properly accounted for: L p k oc 
*pk °' 5 j where L p k and £ p k are the flare peak luminosity 

6 In the time interval 60 s < t < 2 X 10 5 s the 90% contour of 
Fig- El panel (a) is well approximated by: 5.6 X 10 48 [(2^o) _51 + 

(25o) _1 ' 3 ] ' 5 [ 1 + (3i^2) 2 ' 2 ]~ ' 3 - For P anel ( b )' the approximation 
reads: 4.5 X 10" 3 [(|)- 5 - 5 + ( *)-2-6]0.3 
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Figure 4. Average flare luminosity curve for the enlarged sample 
of 55 flares belonging to 29 different GRBs with known red-shift 
calculated following the prescriptions by Lazzati et al., 2008. la 
uncertainties affecting the average value have been shown. Red- 
dashed line: best-fit power-law with index a = —2.6 ± 0.1. As in 
Lazzati et al. 2008 the first data point has been excluded from the 
fit. Blue-dotted line: reference a = —1.5 power-law; the normali- 
sation has been chosen so as to intercept the second point shown. 
Corresponding accretion-rate units are provided on the right y- 

axis: Liso^bcam = f?ac C A/c 2 , with Slbcam = 0.01, JJacc = 0.001. 

The orange hatched area approximately marks the boundary be- 
tween the NDAF (neutrino-dominated accretion flow) and ADAF 
(advection-dominatcd accretion flow) regimes when the neutrino 
cooling switches off (see e.g. Narayan et al. 2001; Chen & Be- 
loborodov 2007). 

and flare peak time as obtained from a lNorris et al.l (|200ot ) 
profile fit. 

In the following we therefore focus our attention on the 
30 s < t < 1000 s interval of time, where the flare luminos- 
ity function (L) oc t~ 2 ' 7 . In Sect. 16.11 we first investigate 
the source of the discrepancy between our results and those 
reported in L08. Possible physical interpretations in the con- 
text of both accretion and magnetar models of the central 
GRB engine are discussed in Sect. I6.ll and l6.2l 

6.1 The L oc t~ 2 7 regime 

The GRB average X-ray flare luminosity function was first 
analysed by L08 starting from a sample of 24 flares coming 
from 9 long GRBs and 1 short GRB. The mean luminos- 
ity was found to decline with a much shallower power-law 
in time: (L)los oc £ _1 - 5±0 ' 2 . This scaling was obtained av- 
eraging the contribution of single flares of assumed square 
shape over time-scales longer than the flare duration. We 
investigate the source of the discrepancy ((L)los oc i" 1 ' 5 vs. 
(L) oc t~ 2 ' 7 of our work) below. We follow the prescriptions 
of L08 and apply their method to the sample of early time 
flares of C10 supplemented by late time flares with known 
redshift: the enlarged sample includes 55 flares detected in 
29 GRBs. We do not assume a flare sq uare shape and used 
instead the best fitting parameters of a iNorris et all (|2005n 
profile: however, the choice of a particular flare profile is 
unlikely to bias the final result, as already noted by L08. 
The interval of times of integration of the mean luminos- 
ity function have been estimated from their Fig. 1. Fig- 
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ure [4] shows the final result: the average flare luminosity 
function is best modelled by a power-law with best fitting 
index a = — 2.6 ± 0.1. The shallower i -1 ' 5 scaling is also 
represented for comparison. Errors are computed following 
the standard theory of error propagation. This result con- 
firms the findings of Sect. with a completely independent 
method. The reason for the discrepancy lies elsewhere. 

We further check the consistency between the L08 
method and ours, measuring the flaring activity evolution 
of GRB 050904. In this case we find L oc t-°- 95 ± 05 which 
is consistent with L oc ^- 10 ± - 3 quoted by L08. The flaring 
evolution of GRB 050904 is much flatter than the average 
behaviour oc t~ 2 ' 7 . Prompted by this result we investigate 
the temporal evolution of the flaring component in the sub- 
sample of GRBs of Table Q] showing multiple episodes of 
flaring emission. The subsample contains 10 elements (black 
dots in Fig. [5]): the mean flare function decay index is —1.8. 
This suggests that multiple flare GRBs have a flatter than 
average flare luminosity function. Interestingly, this subsam- 
ple also shows a flatter than average steep decay: the me- 
dian temporal decay index is found to be evsteep = 1.7 to be 
compared to a s tee P = 2.5 of the remaining 34 GRBs of the 
original sample. This topic is further explored in Sec. 16.21 
The average number of flares per GRB in the L08 sample is 
(nf ) = 2.4 while our enlarged sample contains a wider pop- 
ulation of single-flare GRBs: (n/) = 1.9. In particular, 60% 
of the GRBs of L08 contains more than one flare, while only 
25% of the enlarged sample shows more than one episode 
of activity. The discussion above leads to the conclusion 
that the choice of a sample biassed towards multiple-flare 
GRBs (which was the wider sample of flares with known 
redshift available at the time of writing), led L08 to deter- 
mine a flatter mean luminosity temporal scaling. In agree- 
ment with this picture, in the two GRBs with multiple flares 
for which L08 determined the flare function evolution they 
obtained a = -1.0 ± 0.4 (GRB 050803) and a = -1.0 ± 0.3 
(GRB 050904): both values are flatter -even if consistent- 
than the average —1.5 ± 0.2 decay index derived from their 
entire sample. 

While we benefit from a higher statistics which makes 
our sample more representative of the entire bursts popula- 
tion, the same choice of GRBs with easily recognisable flares 
in their afterglows is possibly introducing a bias which is dif- 
ficult to quantify. The analysis of the flare avoidance region 
of Sect. [SI which properly accounts for the fraction of un- 
detectable flares is the safest approach in this respect: this 
analysis indicates the average flare luminosity decay to be 
steeper than oc t" 1 ' 8 at the 90% c.l. for t < 400 s. 

We finally investigate if the different prompt luminos- 
ity (and consequently afterglow luminosity) of the 44 GRBs 
of Table Q3 could bias the flare (L) temporal scaling. Re- 
normalising the flaring component of each GRB by the av- 
erage prompt luminosity during the burst Tgo duratiorQ the 
(L/Lpr 0m pt) is best fitted by a power-law decay with index 
a ~ -2.5 for t < 1000 s. 

The oc t~~ 2 ' 7 behaviour of the average flare luminosity 
has been obtained in 4 completely independent ways: 

(i) Equation [T] applied to 44 GRBs; 

7 Note that this is different from what displayed in Fig. [2] where 
the x-axis is in T90 units. 



(ii) Equation [T] with re-normalised flaring contribution; 

(iii) I/ P k vs. t p k relation from C10; 

(iv) Application of the L08 method to the enlarged sam- 
ple of 29 long GR Bs with known redshift fitted with a 
iNorris et all (120051 ) profile. 

It is important to underline that the oc t~ 2 ' 7 behaviour rep- 
resents the average scaling of the GRB flaring activity: the 
decay of the flaring component in a particular GRB can con- 
siderably differ from the average relation, as shown by Fig. [5] 
(see also Fig. [2l upper panel, where the yellow hatched area 
marks the 90% sample dispersion). Any model aiming at 
explaining the flare phenomenology is required to provide 
an explanation for the observed scatter, as well. Different 
physical interpretations are discussed below. 



6.1.1 Accretion-powered flares scenario 

In the case of a thin disljf] with constant viscosity for 
t ^> t a „ (where t a „ is the viscous time-scale) the accre- 
tion rate approaches the as ymptotic regime m oc t~ 1,25 
jFranck. King fc Raind Il 985l A similar slope is obtained 
by Cannizzo et al.l |l99fj) performing numerical simulations 
using the a-viscosity pres cription of [Shakura fc Svunvaevl 
1 19751 : rh oc t~ 12 . Recently iMetzger et all |2008r ) computed 
the viscous evolution of an efficiently v— cooled, isolated 



ring of materia l unde r similar assumptions as the study of 
ICannizzo et alj (|l990T l: they report a self-similar behaviour 
rh oc t -4,/3 . The t~ 2 ' 7±0 1 average flare luminosity power-law 
scaling obtained in this work is much steeper than the theo- 
retical expectations listed above: an ad hoc non-linear rela- 
tion between the luminosity and the outflow accretion rate, 
with the efficiency of conversion of mass accretion into lu- 
minosity decreasing as L/rh oc i -1 ' 5 , is re quired to p r eserve 
the thin disk draining accretion scenario . iLee et all (|2009T ) 
demonstrated that the presence of powerful winds driven by 
the recombination of a— particles into nucleons and launched 
from the disk surface could lead to a substantial deviation 
from the t~ A ^ regime, with rh entering an exponential-like 
decay. This regime is unlikely to extend up to t ~ 100— 1000 
s fsee lLee et alj|2009t their Fig. 1). Alternatively, a tempo- 
ral evolution of the flare opening angle causing later flares to 
be less beamed than earlier ones would restore the ~ —1.5 
slop^f). However, most of the beaming evol ution has been 
prove n to take place at early times (see e.g. iMorsonv et all 
120071 and references therein) and is consequently unlikely to 
cause an overestimation of the flare luminosity steepness. 
We therefore consider alternative scenarios. 

Fallback of the stellar envelope that did not reach the 
escape velocity during the explosio n can provide th e source 
of accreting material at late times: IChevaherj [jl989h showed 
that in this case the accretion rate is expected to decline 
as rh oc t~ 5 / 3 which is too shallow to explain the ob- 



served t 



flar e luminosity scaling. Th e process was 



later investigated by Ma cFadven et all (|200lf ): these authors 
simulated the fallback of the stellar envelope following the 



8 Note however that at the high mass accretion rate expected in 
the collapsar model a geo metrically thick disk is expected (see 
e.g. Kumar et al. 2008a.b: iNaravan et al1l200ll) . 

9 In the present work each plot refers to the isotropic equivalent 
luminosity of the X-ray flares at any time t. 
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failure of the Shockwave to unbind the star and recovered 
the rfiRmin oc t~ 5 ' 3 dependency (where mm* is the ra- 
dial fallback rate through their inner numerical boundary 
Rmin = 10 9 cm). However, departures from the t~ 5 ^ 3 scal- 
ing are expected if the infalling material passes through an 
accretion shock developing at R < 10 9 cm or/and if the mass 
fall-back rate decreases suddenly at some time t. These pos- 
sibilities are discussed below. 

Kumar et al., (2008a,b) suggested that the prompt and 
X-ray afterglow reflect the modulation in the rate of cen- 
tral accretion of a rotating progenitor with a core-envelope 
structure onto a black hole. A thick disk configuration is 
assumed. According to their analytical treatment the con- 
dition m oc t~ 3 can be basically achieved in two ways: 
first, in the case of rapid accretion, (£ acc C t), m is ex- 
pected to track the temporal evolution of the mass fall-back 
rate mfp . Assuming a 14 G RB progenitor star (model 
16T1 of IWooslev fc Hegerl 120061 ) . Kumar et al., (2008a,b) 
expect riifb oc t~ 3 (which directly translate into m oc t~ 3 ) 
when the last ~ 0.5 Mq of the star near the surface is ac- 
creted. A second possibility arises when the stellar collapse 
leaves behind an accreting disk and no further mass is being 
added: rhfb = 0. The black hole accretion rate i s found to 
scale as: m w m i~ 4(s+1)/3 jKumar et al.ll2008bl . their Eq. 
43j3- The steepest decline allowed by these conditions is 
m oc t -8 ' 3 , suggestively close to the average flare luminos- 
ity t~ 2 ' 7 temporal decay. This regime establishes for s = 1, 
when a connection-d ominated accretion flow (CDAF, see e.g. 
iNaravan et al.| [2001 and references therein) sets in and the 
accretion rate decreases linearly with disk radius. In reality 
both the steeply falling density profile of the GRB progeni- 
tor and the loss of accreting matter in a wind (which is the 
signature of a CDAF or ADAF scenarios), are likely to play 
a major role in decreasing the fraction of mass which effec- 
tively makes it to the black hole. These two effect coupled 
together are able to account for Lj e t oc t~ 3 or steeper, for 
s > 0.5 (see lKumar et al.ll2008bl . their Fig. 4). 

This scenario w as further explored and confirmed by 
iLindner et al.l d2010h. Using axisymm etric hydrodynamical 
2D simulations. ILindner et all (|2010T ) showed that a steady 
accretion rate rh ~ 0.1 — 0.2 Mqs -1 giving rise to the 
prompt 7-ray phase is followed at t > 40 s by a sudden and 
rapid power-law decline of the central accretion rate which 
is responsible for the prompt-to-steep decay X-ray afterglow 
transition. The steep decline is triggered by a rapid outward 
expansion of an accretion shock through the material feeding 
a convective thick accretion disk and lasts ~ 500 s. During 
this phase the central engine is supposed to be active and the 
accre tion rate is found to decay as m oc t~ 2 ' 7 (|Lindner et al.l 
Again, we find the average flare luminosity function 



to show a very similar scaling (oc t 



for t < 1000 s) 



implying a linear relation between the luminosity and accre- 
tion rate: L ~ ri^cm. The efficiency ?7acc of mass rate to jet 
luminosity conversion depends on many parameters; among 
these, the black hole spin a is likely to play a major role. 
In Fig. U we use ?7acc = 10 -3 which corresponds to a = 0.75 



10 This result is valid for t 2> iacc(fn) where i aC c is the viscous 
time scale of accretion. The scaling of the accretion rate with ra- 
dius has been parametrised as m oc (r/r^) s being the accreting 
disk radius. 
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Figure 5. Best fit power-law index of the temporal decay of the 
X-ray flaring activity vs. continuum (i.e. steep decay) power-law 
index for the subsample of GRBs listed in Table[T]with extended 
flaring emission (black dots). GRBs which show multiple X-ray 
flares but otherwise lack of the redshift measurement are indicated 
with red stars. These events are not part of the main sample of 
Table □ 



following the prescriptions of McKinncv (2005). The main 
limitation of this study lies in that ILindner et al.l (|2010h do 
not simulate the fluid magnetic field which could alter the 
m(t) dependence. 

6.1.2 Magnetar models 

Pulsars arise from core-collapse. A subset of them are 
inferred to host magnetic fields B ~ 10 14 - 10 15 G. 
The rotational energy of a rapidly rotating magnetar 
can be extracted to p ower a strong relativistic jet (e.g. 
iBucciantini et all l2009h that would give rise to the prompt 
emission. The energy release in radiation of this system at 
late times (100 s < t < 1000 s) is subject to great uncertain- 
ties (T. Thompson, private communication). However, the 
output of radiation is likely to be between two limiting cases: 
dipole radiation and monopole radiation. In the former case 
L oc (1+t/T)" 2 which implies L oc t~ 2 for t > T, where T is 
the c haracteristic spin-down time-scale (e.g . Wheeler et al.l 
2000). In the latter, L oc e 21 ^Thompson et alj2004|) . Based 
on the braking indices of observed pulsars in the Galaxy, 
L oc t~ 2 ' 3 is expected (|Bucciantini et al1l2006t ). Addition- 
ally, the interaction with the stellar material can be respon- 
sible for shallower slopes detected in multiple-flare GRBs 
l|Lazzati et al.ll2008h . 



6.2 The variability issue 

Flares are observed as episodic, large-scale amplitude vari- 
ations in the light-cu rves with a typical At/t ~ 0.1 
l|Chincarini et aUl2010l ): this automatically raises the ques- 
tion of the physical source of the flare variability. The pres- 
ence of a correlation between the properties of the prompt 
and flare emission within the same GRB could in principle 
shed light on the mechanism powering the flaring activity. 
However we found no correlation between the prompt lu- 
minosity and the total energy and luminosity of the flaring 
component; the same is true if one were to consider the 15- 
150 keV prompt energy; no correlation has been found be- 
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Figure 6. Upper panel: average flare luminosity function (black 
dots) compared to the average luminosity of the continuum 
underlying the flare emission (blue stars) in the time interval 
40s < t < 1000s. The best fit of both components is also shown. 
Inset: flare-to-continuum flux ratio for the sample of early time 
flares of C10. The dashed line marks the median value of the dis- 
tribution, which is AF/F = AA. Lower panel: flare-to-continuum 
ratio as a function of time. The median value of the distribution 

IS (iflare)/(icont) = 4.7. 



tween the prompt luminosity (or energy) and the peak time 
of the last flare; the peak time of the last flare is also not cor- 
related with the prompt Tgo ; the number of prompt pulses 
is not the key fac tor determining the nu mber of flares to ap- 
pear at t > 50 s (|Chincarini et al.ll2007l ); finally, a hint for a 
prompt-fluence vs. fiare-fiuence correlation was recently re- 
ported by C10: however, the weak correlation is mostly due 
to the presence of two short GRBs. The conclusion is that, 
while flares do share with prompt pulses several key obser- 
vational properties (a notable example is the lag-luminosity 
relation of M10), at the moment it is not possible to infer the 
flaring activity of a burst starting from its prompt emission. 

Flares are under-represented in simple power law X-ray 
afterglows (M10). The X-ray afterglow morphology - flaring 
activity connection is here further explored starting from 
the findings of Sec. 16. II multiple flare GRBs have on average 
flatter flare luminosity functions (i.e. oc t~ a with a < 2.7). 
This, together with the observation that the flare detection 
threshold of Fig. [21 upper panel, closely tracks the temporal 
decay of the average flaring activity for t < 300 s, brought us 
to consider the possibility of a correlation between the X-ray 
flare luminosity decay and the underlying continuum emis- 
sion decay. The average continuum is shown in Fig. [6] the 
best fitting power-law index cY s tee P = 2.8 ± 0.1 is consistent 
with the flare function decay oc t -2 Suggestively, this 

happens for t < 1000 s: at later times the continuum is likely 
to be dominated by the shallow decay component instead of 



the steep decay emission (Figure [6] lower panel, shows a de- 
creasing flare-to-continuum ratio around t ~ 400 s due to the 
progressively increasing contribution of the shallow decay 
component to the continuum). During the first ~ 1000 s the 
typical flare-to-continuum ratio is Lflare/£ s tecp = 4.7 (me- 
dian value) albeit with a large scatter. Using the data from 
C10 we obtain a similar median value -Ffiare/-Fsteep = 4.3 
(Fig. [6] upper panel, inset) . The physical mechanism pow- 
ering each flare emission is therefore required to release an 
average amount of energy -Efl are : 



-Eflare(£) ~ Z/fl are Ai « tL B t ee p(t) 



(5) 



where At is the flare duration. The median value At/t = 
0.23 from C10 has been used. 

The average continuum vs. flare function temporal be- 
haviour of Fig. [6] results from the presence of a correlation 
linking the steep decay flux evolution to the flaring activ- 
ity within individual GRBs. Modelling the steep decay and 
the flaring activity of each burst with decaying power-laws 
of index a s t CC p and anare, respectively, we obtain the result 
drawn in Fig. eO 

flares seems to be linked to the contem- 
poraneous steep decay flux evolution in a way that causes 
flatter flare luminosity functions to be associated to more 
gradual steep decays. 

Generically speaking, these observations are consistent 
with a model where a first physical mechanism is responsible 
for the steep decay continuum (without flares), while mech- 
anism 2 powers X-ray flares. The presence of the a s tcc P vs. 
Qflarc relation suggests that the two mechanisms are in some 
way related and that the sporadic appearance of mechanism 
2 is triggered by some properties of mechanism 1. Insta- 
bilities affecting mechanism 1 can in principle provide the 
source of episodic releases of energy manifesting as flares. 
If this is the case, instabilities are likely to be triggered by 
some physical quantity related to the decay of continuum 
flux at time t: this would explain why late-time flares are 
so rare (only ~ 5% GRBs show clear flaring activity for 
t > 1000 s), but also the paucity of flares in GRBs with sim- 
ple power-law decay X-ray afterglows. In those GRBs, the 
steep decay is likely to be hidden by a conte mporaneous but 
physic ally different emission component (see lMargutti et al.l 
l2010al for a detailed analysis of the two emission components 
in GRB 081028). We speculate that, given the remarkable 
similarity between X-ray flares and prompt pulses, it is pos- 
sible that both mechanisms also operate during the 7-ray 
prompt emission producing slowly varying (mechanism 1) 
and fast varying (mechanism 2) emission components. The 
presence of both long (several seconds) and short (« ms) 
variability time-scales in the same b urst is a known fea ture 
of t he GRB prompt emission (e.g. iNorris et al.l ll99ol but 
also IVetere et ai]|2006l . iBorgonovo et alj|2007ft 

L08 concluded that the outflow interaction with the 
stellar envelope cannot produce flares out of a continuous 
flow. In the following we analyse various mechanisms able 
to produce variability in the central engine release of en- 



11 Only GRBs with multiple flares can be used. Since the knowl- 
edge of the redshift is not necessary, we supplemented the present 
sample with additional 18 Swift-GTtBs with multiple flares and 
good data coverage (red dots in Fig. [5jl. These bursts are listed 
in Table [ATI 
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ergy in the context of accretion (Sect. [6T2. 1|) and magnetar 
models (Sect. [6X2]). 



6.2.1 Variability in accretion models 

In the context of GRB accretion models, a flare corresponds 
to a sudden increase in the jet luminosity as a consequence 
of an abrupt change of the mass accretion rat J* 2 !. Generally 
speaking, large rh variations arise if the disk develops a ring- 
like structure: flares would be associated to the accretion of 
blobs of material initially located at various radii which sub- 
sequently evolve on the viscous time scale. This class of mod- 
els would naturally account for the observed duration-time 
scale correlation and dur ation- peak luminosity anti c orrela - 
tion (CIO) as shown by iPerna. Armitage fc Zhand (2006). 
Different classes of disk instabilities could lead to this sce- 
narico 

Viscous instability arises for drh/dTi < where S is th e 
disk surface density (see e.g. iFranck. King fc Raine|[l985h . 
When this condition is satisfied more material will be fed in 
those regions of the disk that are more dense while mate- 
rial will be removed from less dense areas, so that the disc 
will likely break up into rings. Viscous and thermal instabil- 
ities (the latter taking place for dQ~ /dT < dQ + /dT being 
Q + and Q~ the heating and cooling rates respectively) have 
been invoked to explain Dwarf Novae outbursts: in this case 
the instabilities give rise to a limit-cycle behaviour (see e.g. 
ICannizzo et al.lll99ct ) and do not lead to a total disk break 
down. While the rh required to explain the Dwarf Novae out- 
bursts luminosity is completely different from the conditions 
expected to hold in GRB hyper-ac creting disks, we note 
that outside-in instability bursts (see lFranck. King fc Rainel 
1 19851 and references therein) qualitatively share some obser- 
vational properties with GRB X-ray flares/prompt pulses: 
they have rapid rise -slower decay profiles and rise first at 
longer wavelength (see M10). The stability conditions for 
GRB hyper-accreting disks ha ve been thorough l y stud ied: 
the steady state disk model of |Pi Matteo et all 1 20021 ) led 
the authors to conclude that the accretion flow is both ther- 
mally and viscously stable for a variety of rh values sup- 
posed to give rise to the GRB phenomenology. However, 
X-ray flares suggest that the GRB engine is long -lived: a 
time-dependent comp utation is therefore r equir ed. The time- 
depen dent studies of Ijaniuk et alj 112004 ) and IJaniuk et al.l 
(|2007h revealed the thermal-viscous instability to be an in- 
trinsic property of the innermost disk radii for torus den- 
sities ~ 10 12 gcm -3 ; the GRB disk is not stabilised but 
rather breaks down into rings leading to several episodes 
of dramatic accretion on the viscous time-scale of each 
ring. However, large hyper-accretion rates are required for 
the instability to set in: rh > IOMqs -1 . This value is 
~ 2 order of magnitude higher than the typical accretion 
rate invoked to explain the prom pt GRB phase whic h is 
^prompt ~ 0.1 — 0.2 Mqs -1 (see iLindner et ail |2010| and 
references therein). 



12 The efficiency r? acc = ?? aC c(a) of the process is unlikely to un- 
dergo abrupt temporal changes (see e.g. Kumar et al. 2008a, b). 

13 The importance of hydrodynamical instabilities in collap- 
sar disks has bee n recently demonstrated by the studies of 
iTaylor et ahl feoich . 



Both the black hole spin and the presence of large- 
structure magnetic fields can deeply modify the accreting 
torus structure: for a rapidly spinning black hole, the disk- 
black hole magnetic coupling transfer rotational energy from 
the black-hole to the inner disk where viscous-thermal insta- 
bility is known to arise for a Schwarzschild black hole. The 
magnetic coupling results in reducing the r h value for which 
thermal and viscous instabilities establish: iLei et ail (|2009h 
showed that the disk becomes unstable in its inner region for 
accret ion rates as low as rh > 0.1 Mq s _1 ; Ijaniuk fc Yuan! 
(2010) found that for a Kerr black hole instability in the 
inner edge may be effective for rh ~ 0.5 Mqs - , the ac- 
tual value depending on the viscosity parameter a. While a 
time-dependent analysis is required also in this case, these 
works prove that viscous and thermal instabilities could in- 
deed arise for rh w rh promp t. 

Alternatively, the ring-like structure of the accretion 
disk could be the outcome of a gravitational instability oc- 
curring in the disk outskirt s. Gravitationa l instability arises 
when Q = ^c./ttGE < 1 l|Toomrelll964j l and is known to 
lead to fragmen tation in the outer reg ions of active galactic 
nuclei (see e.g. Shlosma n et al.l 1 19901) and of young stellar 
objects (see e.g. lAdams fc Linlll99ST ). Hyper-accreting GRB 
disks in the steady state regime have been proven to un- 
dergo gravitational instabilities when rh > 2 Mq s _1 and 
radii r > 30 r g (r g = 2GM/c 2 ) under the assumption of a 
viscous parameter a = 0.1. The threshold to gravitational 

r > 15 r g for a — 0.01 



instability lowers to rh > 0.2 Mq s 

llChen fc Beloborodovl 120071 , but see also |Pi Matteo et all 
I2002D . This mechanism has b een suggested as possible 



source of the GRB X-ray flares by Perna, Armitage & Zhang 
(|2006h (see however IProga fc Zhang|2006l : |Piro fc PfhafeOO? 
for a critical view). 

Whatever the origin of the disk fragmentation^ and/or 
ring-like structure is (i.e. viscous, thermal or gravitational), 
the total energy released by each flare is likely to be propor- 
tional to the total mass of the fragment mi which gives rise 
to the burst of accretion. Writing -Efl are ~ Lfl are At, where 
At is a measure of the flare duration, since Lfi arc oc t~ 2 ' 7 and 
At oc t (C10), the average energy of a flare is found to scale 
as -Efl ar o oc f, -1 ' 7 for 30s < t < 1000s. In particular, using 
the best-fitting average flare luminosity calculated in Sec. [3] 
and the flare width vs. peak time relation (w ~ 0.2 x tpjj 
from C10 we find: 



£ is o, fla rc(t) w 1.3 x 10 54 (t/s)- 1 - 7 (ergs" 1 ) 



(6) 



where Siso, flare refers to the isotropic energy emitted in the 
1 — 10000 keV energy band. Equation[S]is valid in the time in- 
terval 40 s < t < 1000 s. A multiplicative factor of 2 has been 
applied to account for the limited 2.2-14.4 keV X-ray energy 
band in which Lfl aro has been computed 1 "^!. Our observations 
would therefore imply the fragments mass mf(t) oc i?fl are (t) 
to scale as: 



14 Note that in the King et al. 2005 picture is the stellar core 
which undergoes fragmentation giving rise to a two-step collapse. 

15 The exact value of this factor depends on the values of the pho- 
ton indices of the Band function describing the spectral energy 
distribution of flares. For the typical low energy and high en- 
crgy spectral photon indices «b = — 1> /9b = —2.5 (Kane ko et al.1 
120061) . and spectral peak energy -E pca k ~ a few keV, this factor is 
fs 2. 
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m,f(t) oc t 



(7) 



Since the arrival time t of a blob of material initially at a 
distance 7? is likely to scale as a positive power of the ra- 
dius t oc B? with > (0 = 3/2 in the case of advection 
dominated flows. [Perna. Armitage fe Zhang : 2006, their Eq. 
2), Eq. [7] suggests that fragmentation taking place at larger 
radii and/or later times gives rise to lower mass objects. This 
is qua litatively consistent with IPerna. Armitage fc Zhana 
(2006), their Eq. 5: however, the actual mass distribution 
of the fragments strongly depends on the local disk proper- 
ties and can be addressed only via numerical simulations (R. 
Perna, private communication). Multiple flares GRBs would 
correspond to a flatter m/ vs. t (and m/ vs. R) dependence 
as implied by the results reported at the beginning of Sect. 



Yet another possibility giving rise to variability in ac- 
cretion models is a modulation of the mass fall back rate ma, 
as proposed by Kumar et al. 2008a,b: while the underlying 
physical mechanism has still to be understood, this mech- 
anism is potentially able to give rise t o sharp variations of 
m provided that £ acc <C t (as shown bv lKumar et al.l l2008b. 
their Fig. 6 and 7) and naturally accounts for the oc t~ 2 ' 7 
behaviour as discussed in Sect. 16.1.11 A possible realisation 
of this process is given by gravitational fragmen tation of 
the f all back material into several bound objects (|Rosswod 
120071 ). 

Magnetic fields can suppress disk fragmentation 
|Baneriee fc Pudritj hoodD , give rise to magneto-rotational 
instabilities (MRI) and strongly modify the dynamics of ac- 
cretion both in magnetars and black-hole sy stems (see e.g. 
IZhang fc Dalll2009l ;]p" roga fc Begelmanll2003l and references 
therein). In particular, a variable output may result from 
a modulatio n of the accretion ra te by the magnetic-barrier 
and gravity l|Naravan et al.ll2003T l . As mass is being accreted 
on to the black-hole, the magnetic flux is accumulated in 
the inner region, causing the accretion through the torus 
to be repeatedly stopped and then restart ed. When accre- 
tion resumes, an X -ray flare is produc e d (|Proga fc Zhana 
120061 ). According to lSpruit fc Uzdenskvl (|2005l ) the accumu- 
lated magnetic field B can support the gas against gravity 
until the radial magnetic force F m is of the order of a few 
percent the gravitational force F g . The force balance yields 
a minimum B oc rh 1/2 . Since the magnetic flux captured by 
accretion also depends on rh, then it is reasonable to expect 
the systems characterised by a flatter m(t) decay (mani- 
fested as a flatter afterglow steep decay) to be able to meet 
the F g ~ F m requirement more easily and repeatedly (i.e. a 
number of flares are expected). This scenario would there- 
fore naturally account for the afl are vs. a a tee P relation of Fig. 

El 



6.2.2 Variability in magnetic models 

This section concentrate on scenarios wh ere the variabil - 
ity completely depends on magnetic effects. iGianniosl (|2006l ) 
suggested that flares can be powered by magnetic reconnec- 
tion triggered by the deceleration with the external medium. 
However, the afl arc vs. a s tccp relation (Fig.(5|, would require 
the steep decay to be intimately connected to the decelera- 
tion phase of the fireball: Swift observations instead favour 
an "internal" origin of the steep decay, being the steep de- 




12 3 
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Figure 7. Flare peak flux to continuum ratio for the sample of 
early-time flares of C10. The distribution can be fitted by the su- 
perposition of two gaussian profiles with the following best fitting 
parameters: x\ = 0.50 ± 0.03, 01 = 0.52 ± 0.03; x 2 = 2.10 ± 0.08, 
0-2 = 0.36 ±0.09. 



cay either interpreted as the tail of the prompt emission (e.g. 
IWillingale et al.ll2010l for a recent study) or the result of the 
prolonged engine activity (e.g. Kumar 2008a,b). 

Alternatively, in the context of magnetar models dif- 
ferentially rotating millisecond pulsars can provide an en- 
gine able to repeatedly store and release energy. Differen- 
tial rotation causes toroidal magnetic fields to be repeatedly 
wound up to ~ 10 17 G and then pushed to and through the 
pulsar surface by buoyant forces: this allows the neutron 
star s pin energy to be emitted in powerful bursts of pulsa r 
wind (|Kluzniak fc Ruderman|[l998l ; iRuderman et al. 2000), 
an extreme but transient realisation of the Usov |l992h" pul- 
sar. This mechanis m has bee n invok ed to explain the flare 
phenomenology by iDai et al.l |2006T l. The re-windup time, 
i.e., the time between subbursts is anticorrelated with the 
pulsar differential rotation AQ; however, the luminosity of 
the subbursts (manifesting as flares) is expected to roughly 
scale as ~ O(t) 4 , bein g Q the neutron star rotation rate 
l|Ruderman et al.ll2000l . their Eq. 19 and 20). If the pulsar 
rotation is also the ultimate source of energy which pow- 
ers the steep decay, then it is not unreasonable to expect the 
correlation displayed in Fig. [(J] and Fig. 



7 SUMMARY AND CONCLUSIONS 

We analysed the X-ray afterglows of 44 long GRBs observed 
by Swift with the goal to determine the average flaring com- 
ponent in the common rest frame energy band which is 
2.2 — 14.4 keV. No assumption on the flare functional shape 
has been made. Our work highlights the importance of the 
proper consideration of the threshold of detection of flares 
against the contemporaneous continuous X-ray emission. In 
particular we showed that: 

• The best fit 2.2-14.4 keV average flare luminosity curve 



(L) = 10 



(54.5±0.1) 



(t/s 



-2.7±0.1 



(ergs" 



for 30s < t < 



1000 s (Fig. [2]). A similar scaling was obtained by C10 for 
the subsample of 43 early-time flares with redshift: L p k oc 
t~,f' 7±0 ' 5 . The application of the L08 method to this sub- 



' P k 



sample of burst leads to (L) oc t 



(Fig.©. 



12 R. Margutti et al. 



• In the time interval 30 s < t < 1000 s the typical 1-10000 
keV flare isotropic energy is found to scale as: -Eiso, flare ~ 
1.3 x 10 54 (t/s)" 17 (erg). 

• For t > 1000 s threshold effects related to the pres- 
ence of a continuum X-ray emission underlying the erratic 
appearance of flares start to play an important role. The 
resulting (L) oc ^- 12 ±° 1 j s biassed towards the bright end 
of the flare luminosity distribution at these times. The un- 
biassed (L) is likely to be steeper (Fig. [2]). 

• According to the flare avoidance region analysis which 
properly accounts for the fraction of undetectable flares at 
any time t, the power-law decay index of (L) is steeper than 
-1.8 at the 90% c.l. for t < 1000 s (Fig. [3]). 

• GRBs with multiple-flare have a flatter than average 
flare luminosity function: (L) oc t~ a with 0.6 < a < 3. Par- 
enthetically, this is probably the reason why L08 determined 
a flatter average flare luminosity function (L)los oc t -1-6 . 

• The decay of the continuum closely tracks the decay of 
the average flare luminosity function (Fig. [FJ| for t < 1000 
s, suggesting that the two components are deeply related 
to one another. In particular, within individual GRBs, the 
power-law decay index of the steep decay is positively cor- 
related to the power-law decay index of the flaring compo- 
nent (Fig. [5]) . GRB 100212A is a show case in this respect 
|Grupe et al.ll2O10l '). The typical flare to steep-decay lumi- 
nosity ratio is L a ^ c /L st cc P = 4.7. 

• As a result, the typical flare energy at time t < 1000 s 
obeys the following relation: Efnac ~ Laare(i)At m tL s t eep (t) 

These findings suggest a model where the steep decay is 
produced by some form of activity of the internal engine 
which w ould be required to be still alive at those times (see 
however I Genet fe Gr anot 2009 for a complementary view), 
while flares could be powered by instabilities affecting the 
physical source of energy which gives origin to the steep de- 
cay. This would explain the X-ray flares erratic behaviour. 
In this picture, the shallow decay phase would be due to a 
completely distinct component of emission which progres- 
sively hides both the steep decay and the X-ray flares as 
time proceeds. 

The (L) oc £ _2 ' 7±0 1 has been analysed in the context 
of accretion and magnetar models of GRBs. In particular: 

• According to the hyper-accreting black hole scenario, 
the L oc t~ 2 ' 7 scaling can be obtained in the case of rapid 
accretion (t acc <C t) or when the last ~ 0.5Mq of the original 
14Mq progenitor star are accreted (Kumar et al., 2008a,b). 
Alternatively, the steep m oc t~ 2 ' 7 behaviour could be trig- 
gered by a rapid outward expansion of an accretion shock i n 
the material feeding a convective disk (|Lindner et al.l feoiO'). 

• The rotational energy of a rapidly spinning magnetar 
can be extracted to power a jet whose luminosity is likely to 
be between the monopole (L oc e~ 2t ) and dipole (L oc t~ 2 ) 
cases. L oc t~~ 2 ' s is expected, based on the braking indices 
of observed pulsars in the Galaxy. 

In both scenarios the variability, which is the main sig- 
nature of the flaring activity, establishes as a consequence 
of different kinds of instabilities. In the case of accretion 
models, thermal, viscous or gravitational instabilities could 
either lead to disk breakdown or fragmentation. Our analy- 
sis constrains the mass of the accreting material to scale as 
mf(t) oc t~ 1,7 (Eq. 7). However, the presence of magnetic 



fields gives rise to MRI instabilities and strongly modifies 
the dynamics of accretion: the accumulation of magnetic flux 
during the accre tion can repeatedly st op and restart the ac- 
cretion process |Proga fc Zhangl2006l ). This would account 
for the erratic flare emission while explaining the aflarc vs. 
Osteep relation. Alternatively, differentially rotating millisec- 
ond pulsars provide a viable mechanism where the existence 
of the afl arc vs. a s teep relation can be reasonably explained. 
We note that if the flare origin is linked to the magnetic en- 
ergy dissipation , the flare emission is likely to be polarised 
jFan et alj|2005t ). while a disk fragmentation origin is likely 
to be accompanied by detectable gravitational wave signal 
|Piro fc Pfhall 120071 ). Both signals will be detectable in the 
near future. 

Whatever the mechanism powering the X-ray flare emis- 
sion is, it is extremely difficult to account for the late-time 
(t > 1000 s) flare activity displayed by some bursts using the 
(L) oc t~ 2 ' 7 component: exceptional circumstances leading 
to the revival of the instabilities would be needed. An inter- 
esting possibility is offered by Fig. [7] while the lower edge 
of the distribution of the flux ratio is probably incomplete, 
this figure suggests the existence of two populations of X- 
ray flares (see also C10, their figure 13). A first population 
with flux contrast ~ 5 (which is the one responsible for the 
average flare and continuum behaviour of Fig. [5} and a sec- 
ond population of bright flares with a typical AF/F ~ 100 
but extending up to AF/F ~ 1000. Fig. [6] directly links the 
X-ray flares to the underlying steep decay flux. It is there- 
fore possible that at late times only the small fraction of 
flares belonging to the second population are able to over- 
shine the contemporaneous shallow decay component. This 
would explain why late time flares are so rare. The detailed 
characterisation of the two flare populations is beyond the 
scope of the present work and is left for a future study. 
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Table Al. GRBs with multiple episodes of flaring activity (Fig. [5}. 
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